Estimation and Control of the Operate 
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Part I — Theory 
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(Manuscript received August 31, 1953J 

The dynamic equations applying to the operation and release of electro- 
magnets are derived in a form in which the armature position and the flux 
linkages of the coil are the variables to be determined as functions of time. 
The effective magnetomotive force and pull are taken as given by the static 
magnetization relations. This formulation is valid if the dynamic field 
pattern is substantially that of the static field, and it is shown that this con- 
dition is satisfi.ed if the effective conductance of the eddy current paths is 
small compared with that of the coil or other linking circuits. This is usually 
the case in operation, but is not the case in normal release. 

Approximate solutions are obtained for both fast and slow operation, 
giving the time as related to the mechanical work done, the mass and travel 
of the moving parts, and the steady state power input to the coil. Expres- 
sions are derived for optimum coil design and pole face area, and design 
requirements for fast operation are discussed. 

An approximate expres-non is derived for the release waiting time, the 
initial period of field decay prior to armature motion. The effect on this 
decay of contact protection in the coil circuit is discussed. The slow release 
case is treated in another article in this issue. An analysis is given of the 
armature motion in release, and the results of an analog computer solu- 
tion to this problem reported. These show the effect of design parameters, 
particularly the armature mass, on the velocity attained in release motion, 
and hence on the amplitude of armature rebound. 

INTRODUCTION 

The operate and release timeH of most telephone relays lie in the 
range from 1 to 100 milliseoH. (0.001 to 0.100 sees.). To use these relays 
in telephone switching, in\'o]\'ins eomplex patterns of sequential switch 
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and relay operation, it is necessary: (1) to be able to estimate the operate 
and release times of any relay for which this information is needed in 
circuit studies, and (2) to develop relays capable of meeting specific 
timing requirements, both fast and slow. 

The essential basis for such design and estimation is a dynamic theory 
of the operation of electromagnets, of sufficient accuracy for engineering 
use. The theory presented here is applicable to electromagnets in gen- 
eral, although the applications discussed are those relating to relays. 
The theory presen'ted is approximate, partly because of the difficulty of 
providing a more exact treatment, and partly because simplicity and 
generality ai'e more important for engineering purposes than accuracy 
in estimation, which is in any case subject to correction by measured 
results. The theoretical relations are used alone only in preliminary 
estimation and in the initial stages of development, as discussed in Part 
I of this article. In advanced development and in the modification and 
application of existing structures, as discussed in Part II, the theory is 
used as a guide in the correlation and extrapolation of observed per- 
formance. 

The dynamics of electromagnets involve the concurrent and inter- 
related phenomena of field development and armature motion, and are 
accordingly governed by the two differential force equations respectively 
applying, each containing a coiipling term expressing their reaction on 
each other. These basic equations are formally identical with those of 
electromechanical transducers, such as loud-speakers, but the treatment 
has fittle else in common. The operation of an electromagnet is the tran- 
sient change from one state of equilibrium to another, as distinguished 
from the sustained low amplitude oscillations of the transducer case. 

The coupling terms represent the effect of field energy changes on the 
coil voltage and on the force causing armature motion respectively. In 
the steady state, the field energy is a function of magnetomotive force 
and of armature position alone. In a transient state, eddy currents in 
the magnetic members affect both field energy and the pattern of the 
field. In developing an approximate theory, it is assumed that these 
effects are confined to the total eff'ective mmf, and that the pattern of 
the field is the same as that of the static field: i.e., that the field energy 
associated with any portion of the structure, such as an air gap, is fixed 
by the flux linkages of the coil and by the armature position. The limita- 
tions on the analysis imposed by this simpHfying assumption are dis- 
cussed at the end of the next section. 
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1 The Dynamic Equations 
The notatiou of this article conforms to the list given on page 257. 

THE electrical EQUATION 

The voltage equation for a coil of iV turns linking a field of strength tp 
may be written in the form: 

a-/)/? + A^^ = o, (1) 

where i is the instantaneous current, 7? is the circuit resistance, and IR 
is the constant voltage applied. Writing 5, for the instantaneous mmf 
4ff Ni, and 3^,>j, for the steady state mmf 4Tr NI, the equation becomes: 

di 

where Gi = N /R, the coil constant, or equivalent single turn conduc- 
tance. If the field ip is linked l>y a number of circuits, a similar voltage 
equation applies to each such circuit. By addition of these expressions, 
there is obtained: 

ff - i?. + 47rG ^ = 0, 
at 

where 3^ = ^ffv , the total effective mmf, fFs = ^fF^i , the total applied 
mmf, and G = X!^* > ^^^ total equivalent single turn conductance. If 
the dynamic field has the same pattern as the static field, the instan- 
taneous mmf 5 must equal (Rip, where 6i is the reluctance ^/tp observed in 
static magnetization measurements. The preceding equation may there- 
fore be written in the form: 

6lip - S, + 47rG ^ = 0. (2) 

at 

This relation controls the time rate of development of the flux ip. As 
the coil is the only circuit linkmg the field that has an externally applied 
voltage, Js is simply the steady state applied mmf -iir NI. The conduc- 
tance G includes not only the coil conductance Gc , or N'^/R, but terms 
for all conductive paths linking the field, including the eddy current 
paths. The effective conductance of the eddy current paths is denoted 
Gb . When a short circuited ending or sleeve is used, as in slow release 
relays, its conductance is denoted Gs . Thus in most cases of interest 
G = Gc -\- Ge ; when a sleeve is used, G = Gc -\- Ge -\- Gs . 

The static magnetization relations between £F and (p vary with the 
position of the armature, which is defined by its displacement x from the 
operated position. Hence the reluctance fll is, in general, a function of 
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both X and tp. Thus equation (2) is a differential equation in which tp, x, 
and dip/dt are the variables. With the armature at rest, x is constant, 
and (2) may be solved to give the initial field development in operate, 
or decay in release. When the armature is moving, the variation of <p and 
X is governed jointly by (2) and the mechanical force equation. 

THE MECHANICAL EQUATION 

The forces acting to accelerate the armature and associated moving 
parts, of effective mass m, are the magnetic pull F and the forces exerted 
by the contact and other springs. The total spring force may be written 
as dV/dx, where V is the potential energy stored in the springs, expressed 
as a function of x. The force equation may be written in the form : 

As X measures the gap opening, the velocity dx/dt is positive as the 
gap opens, and negative as it closes. The pull is directed toward the 
closed position and therefore tends to algebraically decrease the velocity. 
The spring force tends to open the gap and to algebraically increase the 
velocity, as V decreases with x, and dV/dx is negative. 

On the assumption that the dynamic and static field patterns are the 
same, the pull F can be evaluated from the static magnetization rela- 
tions, i^ is a function of (p and x, the variables of the flux development 
equation (2). Thus the dynamic performance ia governed jointly by (2) 
and (3), to which the concurrent variation in <p and x with time must 
conform. To complete the formulation there are required explicit ex- 
pressions for (R and F in terms of (p and x. 

HELUCTANCE AND PULL 

As 8ho\vn in a companion article,^ the magnetization relations for 
most electromagnets conform to the simple magnetic circuit of Fig. 1. 
Through the region of linear magnetization, in which the flux density is 
below that producmg incipient saturation, the reluctances (Rn and CH/, 
and the area A are constants. If the magnetization relations conform 
to this schematic, the total reluctance (R is given by: 



(Rl ((Ro -I- 2) 



(R = ^ ^ . (4) 

The constants (Ro , (Rl , and A of Fig. 1 and equation (4) are called 
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the equivalent values of the closed gap reluctance, the leakage rehictauce, 
and the pole face area respectively. Procedures for estimating these 
quantities from the magnet's dimensions and material constants are 
described in the article cited above, while methods for their experi- 
mental evaluation are described in another article in this issue.' When 
equation (4) applies, the pull F is given by: 



F = 



i&L'py 



SttA ( (Rl + Cfio + J 



(5) 



In normal operation, the flux level lies in the region of linear magnetiza- 
tion, and equations (4) and (5) apply through the greater part of the 
period of flux development determining the operate time. In what fol- 




Fig. 1 — Equivalent magnetic circuit of an electromagnet. 

lows, some simplification of notation is obtained by WTJting .Tq for ASio , 
and C/, for ((Ro + (Rl)/(Ro , so that these equations become: 

{Cl- l)(.i-o + .r) 



(i{ ^ 6{ 



F - 



ClXo + X 
(C, - l).ro 

ClXq + X 



8irA' 



(4A) 
(5A) 



Thus the dynamic relations are given by (2) and (3), in which Oi and 
F arc given by (4) and (5) respectively. The magnetic circuit constants 
can be evahiated by the methods cited above, and the other constant 
terms are known or given quantities, with the exception of the eddy 
current conductance Ge . Evaluation of this term requires determination 
of the conditions under which such a constant term can adequately 
represent the effects of eddy currents. 



EDDY CURRENT CONDUCTANCE 



Fig. 2 shows a simple electromagnet with a cylindrical core of length t 
and diameter D. If (it is the rehictance, the magnetomotive force fFc of 
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the winding current results in a flux <pc = fFc/01, uniformly distributed 
across tiie core cross section. Tlie leakage paths cause some longitudinal 
variation in ^, which can be neglected for the present purpose. If the 
flux is varying, eddy currents flow in circular paths with a current den- 
sity j, which varies with the radius r. These give an increment to the 
total mmf varying from zero at the surface to a maximum at the center, 
producing a correspondmg variation in the density of the total flux. 
If '3c is large compared with the maximum magnetomotive force pro- 
duced by the eddy currents, however, the density is nearly constant, 
and its rate of change is nearly constant throughout the core cross sec- 
tion. To the extent that this condition is satisfied, the effect of the eddy 
currents can be determined as follows. 




Fig, 2 — Eddy current paths in the core of an electromagnet. 

The current in a cylindrical shell of radius r and differential thickness 
dr is jl dr. This links a part of the held proportional to the area enclosed, 
or ^r^tp/D^, which varies, by hypothesis, at the same proportional rate 
as the total field. The resistance of the shell is 2Trrp/{t dr.), where p is 
the resistivity of the material. The voltage equation for the shell cir- 
cuit is therefore : 

dtp 



2jrrp 4r »^ „ 



and hence: 



3 = 



2r 
pirD' 



dip 
dt 



The magnetomotive force of the shell is its current multiplied by 4jr. 
This produces a flux increment d^p in the area enclosed inversely propor- 
tional to the reluctance of the tubes of induction mthin this area. This 
reluctance may be taken as inversely proportional to the area, as would 
be the case in a closed magnetic circuit of uniform cross section. Then 
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the flux increment produced by each shell is given by; 

On substituting the preceding expression for the current density j, 

rDt2 

and evaluating the integral / dip, the total flux <ps produced by the 

Jo 

eddy currents is found to be: - 

This is identical in form \\'ith the expression for one of the linking 
circuits assumed in deriving equation (2), with (R^s equal to Se , and 
^.tGe given by t/{2p). Thus to the extent that the assumption of a uni- 
form flux distribution is valid, Ge is given by: 

This expression applies to a round core. A parallel approximation may 
be obtained for a rectangular rod by considering it as made up of rec- 
tangular shells of differential thickness, having their sides in the same 
ratio as the rod. By treating the perimeter of each shell as corresponding 
to the circumference of the circular shell in the cylinder case, with the 
area enclosed corresponding to that enclosed by the circular shell, the 
following expression is obtained for the conductance of a rectangular 
rod whose sides are in the ratio k: 

This shows that the effect of using a rectangular section is equivalent 
to increasing the resistivity of the material. A similar effect is obtained 
by subdividing the core section, as in laminated construction. If the 
core, for example, were made of two similar round rods, equation (6) 
would apply to each, and the effective mmf acting on both would be the 
value of (SiipE given by this expression, with tp in d<p/dt equal to half the 
total flux, so that the resulting expression for Ge would be half that given 
by (6). This argument can be generalized to show that, aside from the 
effect of changes in section shape, the effect of subdivision is to reduce 
Ge in proportion to the number of subdivisions. 

The approximate validity of these expressions for Ge rests on the 
assumption that the eddy current magnetomotive force is minor com- 
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pared with that of the coil or other external circuit. From the derivation 
of (2), the component mmf' s are proportional to the corresponding terms 
in G, so the preceding expressions for Ge are valid only if Ge is small 
compared with Gc(N'/R). In all but exceptional cases, this condition is 
satisfied in operation, and in release with a short circuited \Wnding or a 
slow release sleeve. In these cases, the effect of eddy currents is ade- 
quately represented by a constant Ge term in (2). The exact values of Ge 
applying are in only approximate agreement mth those given by (7) 
and (8), as the derivation of these expressions ignores the variation in 
flux density along the length of the core, and the eddy current effects in 
the armature and return path. 

EDDY CURRENTS IN RELEASE 

In normal release, the winding circuit is open, and the only effective 
magnetomotive force is that of the eddy currents. Evidently, the field 
in the outer layers must collapse almost instantly, while that in the 
center of the core is sustained by eddy currents in paths whose mean 
conductance, per line linked, is higher than that applying to a uniform 

field. 

The relations applying to a closed path of uniform section can be 
formulated in differential form, and the solution for a cylindrical section 
has been given by Wwedensky.^ For decay from an initially uniform field, 
the expression for the flux as a function of time is a series of exponential 
terms with progressively smaller time constants. The first term represents 
the most persistent part of the flux, and represents a field varying from 
zero at the surface to a maximum at the center, comprising 69 per cent 
of the initial field. Its time constant corresponds to an effective value of 
Ge 35 per cent larger than that given by (7). Somewhat similar relations 
must apply to an electromagnet, causing a time variation in the pattern 
of the field, not only radially, as in a (closed uniform path, but in the 
longitudinal variation and in the division of the field between the leak- 
age and armature paths. 

An experimental study of flux development and decay in relays has 
been reported by M. A. Logan.^ His results agree with this discussion in 
showing Ge in (2) to be effectively a constant, provided Ge/Gc is less 
than 0.2, as in the operation of most relays. An empirical expression is 
given for an effective value of Gc + Ge which provides a correction fac- 
tor applicable for small values of Ge/Gc . The results for normal release 
show the field decay to have the general character of Wwedensky's 
solution, and an empirical expression is given which agrees with ob- 
served results. This is primarily of interest in connection with the voltage 
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induced in the coil and imposed on the contact opening tlio coil circuit. 
Because of the changing field pattern, the gap field which determines the 
pull has a dilTerent, lliough similar, decay rate. 

These considerations, as supported by Wwedensky's relations and 
Logan's results, indicate that the rate of pull decay in normal release is 
faster, and only roughly of the same order of magnitude, as that esti- 
mated on the assumption that (2) applies, with values of Oe similar to 
those applying in operation. 

LIMITATIONS OF THE ANALYSIS 

The validity of equations (2) through (5) rests on the assumption that 
the pattern of the dynamic field is essentially that of the static field to 
which the magnetization relations apply. The above discussion of the 
eddy currents indicates that this assumption is valid when Ge is a minor 
term in G. This condition is approximately satisfied in normal operation. 
In open circuit release, however, the condition is not satisfied, and equa- 
tion (2) is only a crude first approximation to the controlling relation. 

The use of (4) as an expression for the reluctance 01 rests on the fur- 
ther assumption that the magnetization is linear, a condition only satis- 
fied in the low density region. The initial and controlling stages of opera- 
tion are usually complete liefore the field passes out of the low density 
region, and (-1) is therefore applicable to the operate case. A diftorent 
expression for (il is required in the release case, as discussed in Section 7. 

2 Character of the Operate Solution 

GR.\PHICAL representation 

Some understaTuling of the relations applying to operation may bo 
obtained from their graphical representation in Figs. 3 and 4. In these 
two figures the path followed by the variables in dyamic operation is 
indicated by the dotted lines, with the dots spaced to indicate equal 
time intervals between 1 hem. Fig. 3 shows the ip versus 5 relation, referred 
to the steady state magnetization curves for various values of x, with 
X = corresponding to the operated position, and x = Xi to the initial 
unoperated position. Fig. 4 shows the dynamic /'' versus .t relation, to- 
gether with the load curve (bounding the cross hatched area V) and the 
steady state ]ni!l rwvvv for the applied mmf Ss ■ 

The flux and pull increase together with the armature at rest at Xi 
until the pull equals the back tension at the pohit 1. In the earlier mo- 
tion, 1-2, the \'elocity is small, and the reluctance (from (4)) changes 
slowly with x so the motion has little effect on the rate of flux develop- 



118 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1954 







.^ 


^•'"^ 


"1 






^ 


0.2 




/ l"""""'''''''''^ 








/ JX» 






/ 




-— - 


0,4 


/ 


/ .^'""^^ 






/ / 


/ *^ 


. 


0.6 


1 / 


/ ^•'"'^^ 


^ 


. 0,8 


// 


//v^^ 


■— - 


5*51; = ^ 


///^ 


1^^^ 






^ 









3^s 



MAGNETOMOTIVE FORCE, J 

Fig. 3 — Field energy relations in the operation of an electromagnet. 

ment. In the later motion, 2-3, the rehuitance changes more rapidly 
with X, and the velocity is high, increasing d(p/dt so as to result in a 
temporary decrease in ^. Operation is complete at 3, with (p and F still 
below their steady state values at 4, which they then approach exponen- 
tially with X = . 

The mechanical work done in operation is represented in Fig. 4 by 
the area under the dynamic pull curve (dotted line), and in Fig. 3 by 
the area bounded by 0-1-2-3-0. This, of course, is less than the work 
that would be done if fF were equal to its steady state value ^s throughout 
the motion, represented by the area under the fF.s curve in Fig. 4, and 
by the loop 0-3-4-5-0 in Fig. 3. In Fig. 4 it can be seen that the work done 
exceeds the static load V by an amount represented by the shaded area 
T, corresponding to the kinetic energy of the armature and the parts 
that move with it. At the end of the stroke this energy is partly dissipated 
in impact, and partly transferred to vibratory motion of the relay and 
its parts. 

The initial flux development, 0-1 is governed wholly by equation (2). 
This same relation dominates in controlling the performance in the early 
travel 1-2, in which the velocity is small, and the reluctance changes 
slowly. In the later travel, where the velocity is high and the reluctance 
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changes rapidly, equation (3) is controlling, mth F corresponding to a 
more nearly constant value of <^. 

METHODS OF SOLUTION 

Formally, the expression for ip given by combining (3) and (5) may be 
substituted in (2), eliminating >p and gi\dng a third order non-linear dif- 
ferential equation in x and (. No analytical solution to this equation is 
known, even for the idealized case of a linear load curve. Specific solutions 
may be obtained by means of computing devices, but there are too many 
parameters involved to prepare a useful library of specific solutions that 
would serve as a general treatment. 

Approximate analytical solutions can be obtained by a step by step 
process, starting from the initial stage of flux development with the 
armature at rest. These solutions are relatively simple in form and ap- 
plication, and provide an anah'tical statement of the relations between 
performance and design parameters, from which conclusions can be 
drawn as to design features and operating conditions favoring a desired 




TRAVEL, 'X — * Xi 

Fig. 4 — Load and jiull relalioiis in the operation of an electromagnet. 
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type of performance. For development purposes, an approximate solu- 
tion indicating such optimum conditions is preferred to a rigorous solu- 
tion in a form too complex for such use. 

Three such approximate solutions for the operate case are described 
below. The first of these, the single stage approximation, is an applica- 
tion of the solution to the flux rise equation (2), neglecting the armature 
motion. It is applicable in cases of slow operation controlled by the rate 
at which the flux development provides sufficient pull to operate the 
load. The other two solutions are more generally applicable, particularly 
for relatively fast operation. One of these, the two-stage approximation, 
is the simpler in form, while the other, the three stage approximation, 
is the more accurate. 

3 Single Stage Approximation 

The initial stage of flux development with the armature at rest ends 
when the pull equals the back tension. During this stage the flux develop- 
ment is governed wholly by (2), with x = xi , the initial gap, and 
Si = (Ri , as given by (4) for x = xi . Writing (Pi(Sii for Js , equation (2) 
becomes : 

4xG [■•' d<p 
CRi Jo Vi — <p' 

which on integration gives : 

4irG 1 

where v = ipl*f>\ , the ratio of the flux attained at time i to the steady state 
flux <^i or fys/(Jli ■ The mmf ratio g^/Js is also equal to v. 

Equation (9) applies rigorously only while the armature is at rest. 
It is, however, also a close approximation for the initial motion, in which 
the velocity is small and x difTers little from x^ . If the rate of flux de- 
velopment is slow ((?/(Ri large) the armature moves slowly mth the 
pull only slightly in excess of the load curve. In this case the inertia 
term m-ctx/<l^ in (3) is minor, and the pull F nearly equals the static 
load dV/dx. Thus the operate time is approximately equal to the time 
required to develop a pull which exceeds the static load at all points in 
the travel. Tliis pull is attained at the just operate current or ampere 
turn value, corresponding to the minimum mmf for static operation, 
ffo - A.ssuming that the armature moves as the flux and pull develop, the 
operate time is that required for 5 to equal Jo , and is given by equation 
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(9) for r = fF()/Js ■ While this is a cmde approximation, it provides 
satisfactory estimates of the time for slow operation. 

EDDY CURRENT CONDUCTANCE DETERMINATION 

Equation (9) is used in one experimental method for the evaluation 
of the eddy current conductance Ge . In this, measurements are made of 
the time at which motion starts for various values of the coil constant Gc ■ 
The latter may be varied by adding series resistance, and the applied 
voltage adjusted to maiiitaiTi the current and hence ^s constant in suc- 
cessive measurements. The time / for motion to start may t>e determined 
from shadowgi'aph measurements. For a constant bac^k tension, the 
pull and hence the value of v or ET/fFs when motion starts is constant in 
these measurements. From (9), / is then directly proportional to Gc + 
Gs , so that a plot of i vs. Gc should be linear, with a negative intercept 
on the Gc axis numerically equal to Gs . Experimentally, an approxi- 
mately linear relation is obtained in this way, provided the values of 
Gc covered are in excess of oGe , in agreement with the discussion of 
Section 1. Values of Gs thus determined are consistent with those obtained 
from other measurements, and in approximate agreement with esti- 
mates obtained from (7) or (8). 

COIL CONSTANT 

In equation (9) Ge and (Ri are determined by the magnetic design. 
The latter, together with the load, determines Jo , the just operate mmf , 
or 4ir(NI)o . As V equals Jo/STs , or {NI)q/{NI), the only quantities not 
fixed by the design and the load requirement are Gc and the steady 
state ampere turn value .¥/. As the square of this latter quantity is 
equal to Gc-I R, equation (9) may be written in the form: 

,.;|5(c,,.+ (i;g»),n-l-. (10) 

(it, \ V-l'R / \ — V 

The steady state power J'R is either determined by circuit require- 
ments, or (chosen with reference to economy of power consumption. 
With I^U fixed, the only independent quantity in (10) is v, which is 
determined by the choice of the coil constant Gc . Neglecting Ge , the 
time is proportional to 

In , 

V- v — V 

which is shown plotted against v in Fig. 5. This figure includes a curve 
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giving corresponding values of In 1/(1 - tO, which may be used in 
specific cases to determine the correction corresponding to the term 

in Ge ■ 

In most cases of slow operation, the operate time is of little im- 
portance, and Gc is chosen to make NI greater than (NI)o in all cases: 
i.e., after allovvdng for possible variations in {A^/)o resulting from varia- 
tions in the load and in the magnetic characteristics. These variations 
result in a variation in v, and the corresponding variation in time is shown 
by the curve of Fig. o. If, however, it is desired to minimize the operate 
time for a given power input, Gc should be choseji to give the value of 
V (0.715) corresponding to the minimum of the curve. As this minimum 
is broad, variations in v, corresponding to those in (A'^/)o , produce little 
change in tne operate time. 

4 Two Stage Approximation 

Unless the rate of flux development is very slow, as assumed in the 
single stage approximation, the pull attained in the early travel is in 
excess of the load, and the kinetic energy T is a considerable part of 
the total work output V + T. As a first approximation to this case, the 
operate time can be computed as though operation occurred in two suc- 
cessive stages : (1) a stage of flux development with the armature at rest 
in the unoperated position {x = Xi), and (2) a stage of motion, in which 
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Fig. 5 — Relations for estimating time of flux development. 
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the flux remains constant at the value attained at the end of the first 
stage. This approximation necessarily over-estimates the operate time, 
as it takes the two stages as occurring in sequence, whereas they ac- 
tually overlap. 

Let /i be the time of the first stage, and /2 the time of the second stage, 
where I2 is the time for the motion from xi to some other armature posi- 
tion .Tn , measuring the completion of operation. The first stage ends when 
the flux has attained some value vipi , where (pi = EFs/cRi , the steady state 
fiux for the open gap, and /' remains to be determined. Then ti is given 
Ity (9) for this value of v. 

By hypothesis, the flux remains constant at ^i^i throughout the second 
stage, and the corresponding pull is then given by (5A) for tp = vtpi . 
The work V -{- T is the integral of F-dx from X] to X2 , given by the 
integral of (5A). Hence: 

"^ Sir A {C^o + x,){Cao -\-X2)' ^ ' 

In this equation Xo/A may he replaced by the expression for (Ro 
given by {4A) for .-r = Xi , and Gc-l'R substituted for (NI)^ in the re- 
sulting equation, which then becomes: 

where Cw is given by : 

— (.To + .nXC^-To + ^Ta) , , 

(C,, - l)Xn{Xi — X2) 

The time t2 of the second stage depends upon the difference between 
the pull and load curves. As a representative condition, this may be 
taken as approximately constant. For such uniform acceleration of the 
effective mass m through the distance Xi — x^ , the kinetic energy T 
equals 2mCxi — X2f/tl, giving an expression for (2 in terms of T. As the 
time ^1 of the first stage is given by (9), the total operate time to , or 
ti + is is given by: 



MGe + Gc) 

in — - 



l_^(2,nfc-^y (14) 



(it, 1 - V 

In this approximate expression for the operate time, v and T are as 
yet undetermined quantities, related by equation (12), in which Gc 
is tlui value of N'/R for the coil used. The two stage o])eration assumed 
in deriving this expression corresponds to the existence of a restraint 
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which holds the armature at rest until it attains the flux level v<pi . The 
approximation will most nearly approach the true relation if the effect 
of this restraint is minimized, and v taken at the value which minimizes 
to . This minimum value can be determined by trial, using the curves of 
Fig. 5 to compute /(, for several values of v, and taking the minimum from 
the resulting plot of lo versus v. 

In development studies, the operate time of interest is usually not 
that for an assigned winding, but that for the winding gi^dng minimum 
time for a given power input. Substituting in (14) the expression for 
Gc given by (12), there is obtained the equation: 

where tg is written for 47rGs/CRi , the eddy current time constant. As 
Gc remains to ho determined, this expression contains two independent 
variables, v and T, which are to be so ehosen as to make ^o a minimum. 
Equating the partial dei'ivative mth respect to 3" to zero, there is ob- 
tained the following ex-pression for the value of T for which the time 
is a minimum: 



T = 2m(;r] - X2) 



2 I v'fR 



^Cw In ^ 



1 -v/ 
On substituting this expression for T, the preceding equation becomes: 

If the value of v which minimizes /o is determined, the resulting value 
of tn is the minimum operate time attainable by optimum coil design. 
In particular, if t^ is negligible, this minimum corresponds, from Fig. 5, 
to V = 0.715, for which 

v^ I — V 

The corresponding value of Gc , the coil constant value making l„ a 
minimum, is given by (12). In this use of equation (12), v and T are 
taken at the optimum values obtained as indicated above. 

OPTIMUM COIL CONSTANT 

As in the case of the single stage approximation, the optimum coil 
constant corresponds to a value for v of 0.715 when Ie is negligible. From 
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the curves of Fig. 5, it can he seen that if ts is not negligiljle, the optimum 
vakiG of V is less than 0.715, and will be smaller the larger the ratio of the 
first term in (15), that in In, to the other two terms. From (12), a reduc- 
tion in V corresponds to an increase in yi, and hence in the coil constant 
Gc (/ R being given). Thus the optimum coil constant for fast operation 
is larger for electromagnets with long cores of low resistivity material, 
for which Ge and (e are large, than for those with short cores of high 
resistivity material, for which these quantities are small. 

The physical significance of the optimum coil constant can be deduced 
from the relations shown graphically in Fig. 4. For a given value of the 
power input PR, an increase in Gc increases both the steady state mmf 
Ja and the time constant of flux development 47rGc/tRi , increasing the 
upper limit to the attainable pull while reducing the rate at which this 
limit is approached. If the coil constant were so small that the area 
under the dynamic pull curve equalled the static work load V, the operate 
time would be infinite. On the other hand, an infinitely large coil con- 
stant would correspond to infinite inductance and an infinite operate 
time. Between these extremes lies the optimum value of the coil con- 
stant, gi\ang minimum operate time. This is a broad optimum, near 
which a change in Ja is compensated by a corresponding change in the 
rate of flux development, so that the realized dynamic pull curve is not 
affected by a small change in Gc . 

In addition to the operate time, the value of Gc affects the final 
velocity of the armature, and thus its kinetic energy in impact with the 
core at the end of the stroke. This energy is dissipated in the relay and 
spring vibration associated with contact chatter. For values of Gc above 
the optimum, the higher inductance reduces the pull in the early travel, 
while the higher value of 5s increases it in the later travel. The net effect 
is to increase both the operate time and the final velocity. Values of Gc 
aljove the optimum are thorefoi'c disadN-antageous not only in slower op- 
eration, but in increased impa(^t energy tending to cause contact chatter. 

FACTORS CONTROLLING SPEED 

Aside from the term in Ie , equation (15) shows that the minimum 
operate time, corresponding to the optimum value of v, is determined by 
the static load V, the inertia load measured by m(xi -~ x^^, the steady 
state power 7-7?, and the constant d, . The latter is given by (13), and 
is the only quantity in (15), aside from (e , which depends upon the mag- 
netic design. In the range of values applying in practice, Cw is deter- 
mined primarily by the leakage factor Cl , measuring the ratio of leak- 
age to useful flux. In most practical cases, X2 is small (zero for complete 
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operation), Cl equals or exceeds 4 ((Jit equals or exceeds 3 (Ro), and aii/aio 
lies in the range from 1 to 3. For 4 < C,. < 10 and 1 < xi/xo < 3, 
1.5 < C'w < 2.0. In practice therefore Cw has a value close to 2.0 for 
most electromagnets. 

The term in V, the second term in (15), varies inversely as the power 
fR, while the third, or inertia term, varies as the cube root of the power. 
Hence the second term tends to dominate at low values of PR, and the 
third term tends to dominate at high values. For a low power input 
therefore the operate time is controlled by the load V, and varies in- 
versely as the power input, while for a high power input, the operate 
tine is controlled by the inertia term, m(xi — x^f, and varies inversely 
as the cube root of the power. In the load controlled case, the time varies 
directly as the load; in the mass controlled case it varies as the cube 
root of the mass and as the two-thirds power of the travel Xi - 3^2 ■ The 
eddy current term is very nearly a constant increment to the other two 
terras, and is therefore relatively more important, the faster the opera- 
tion. 

PRELIMINARY TIMING ESTIMATES 

Within the limits of accuracy to which this two stage approximation 
applies, the effect of the magnet design upon the operate time appears 
only in is and Cw . If the former is neglected, the optimum value of v 
is 0.715, and 

\\n^ =2.5. 

Taking Cw as having the representative value of 2.0, (15) reduces to: 

lOF , / 1 35m(.ri - x^Y ^ . . 

^'^PR'^ \ PR / ■ ^ ^ 

This simple approximation provides rough estimates of the operate 
times attainable with any electromagnet for a given load and power 
input. As an illustration, consider a typical relay spring load involving 
a travel of 40 mil-in (0.1 cm) with a level of spring force of 100 gm (10^ 
dynes), so that V = 10* ergs. Let the effective mass m of the moving parts 
be 10 gm. For a steady state power input of 0.5 watts (5 X lO" ergs sec), 
the two terms of (16) have values of 0.020 sec and 0.014 sec, respectively, 
so that ta equals 0.034 sec. For an input of 5 watts, the two terms become 
0.002 sec and 0.007 sec, respectively, and ta equals 0.009 sec. The neg- 
lected ts term of (15) might amount to an increment of 0.005 sec for a 
solid core of magnetic iron, but would be proportionally smaller for higher 
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resistivity material. In either case, this increment would be of little 
consequence in the low power case, but would materially affect the 
operate time in the high power case. 

5 Three Stage Approximation 

The following analysis provides greater accuracy in estimating operate 
time than the two stage approximation, together with a more accurate 
representation of the relations between performance and the controlling 
variables. It differs from the two stage approximation in treating the 
initial motion as a separate stage of operation. The three stages of 
operation thus become: (I) Increase of the flux to the value at which 
the pull equals the back tension, mth the armature at rest, (2) flux 
development and concurrent armature acceleration, assuming equations 
(2) and (3) to apply, with the approximation that the variations of 
(R and ff with .-c are ignored, and (3) the later motion, which is treated 
in the same manner as the second stage of the two stage approximation. 

Formally, the second stage should be restricted to a very small part 
of the total armature motion, in order to minimize the change in (R 
which is ignored. Relatively minor error, however, is introduced in ig- 
noring the variation in (fi through as much as half the total travel. The 
spring force is taken as constant through the second stage. In the relay 
case, this is approximately true for the initial travel prior to the actuation 
of the contacts, which results in an abnipt increase in the spring load. 
Thus the second stage can be taken to extend to the travel at which 
contact actuation first occurs. If all contacts are actuated near the same 
point in the travel, the end of the second stage coincides with complete 
operation, and the third stage need not be considered. If the contacts are 
spread out, the third stage coincides ^ath the stagger time, or time be- 
tween operation of the first and last contacts. The operate time there- 
fore consists either entirely or principally of the time for the fir.'it two 
stages. 

As before let: 

xi be the initial (open) gap, for which (Jt, = Oti , 

7^1 be the back tension (constant load in the first two stages), 

tc be the initial coil time constant, — ■ =— , 

(Ri R 

47r 
ts be the initial eddy current time constant, — Ge, 

tRi 

Ss be the steady state mmf, ivNI, and <pi the corresponding flux for 

X = a'l , or Ss/is\\ . 
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FIRST STAGE 

The pull for x = xi is given by (oA), which may be written: 

where Ai is given by: 

A, = ( ^f° + "' Y A. (18) 

\Xo{Cl — 1)/ 

Tn the firKt stage the flux increases to vm , where vi is given by (17) 
for F = Fi . Then the time ^i of the first stage is given liy (9) for v = vi, 
or by : 

h = (ic + Ie) In —^ . (19) 

1 — Ul 

The solution for the first stage therefore determines fi and ti . 

SECOND STAGE 

It is convenient to write the equations for the second stage in terms 
of the ratio z defined by: 

t 
tc ~r *E 

where t is now measured from the start of the second stage. Neglecting 
the change in reluctance with armature motion, the ratio <p/ipi or v is 
given by (9), except that for v = Vi at f = the relation must be written: 

y = 1 - (1 - Vi)e-' . (20) 

As <pi = O^a/ffli , or 4TA^//mi , and '^ "" ^ " ^g ' 



2 4x ,2^ 



(f>i = — I R'tc ■ 

Taking the pull in the second stage as given by the expression apply- 
ing for X = xi, and substituting the preceding expressions for v and 
(,31 , there is obtained : 

"- ii^^' -''' -'""''■ <^« 

For the second stage dV/dx is constant at -Fi throughout, and 
(3) reduces to m d'x/df = F - Fi , or: 
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dt^ 2m(RiA: 



((1 - (1 - v,)e-')' - vl). 



Substituting {tc + te) z for t, this equation may be integrated with 
respect to t to obtain the following expressions for the velocity x and 
the travel x in the second stage: 

2m(i{iAi 

The functions /](vi , 2) and/2(yi , s) are the integrals with respect to 2 of 
the bracketed term in the preceding equation. These functions can be 
evaluated from the curves of Fig. 6, which shows them plotted against 
z for several vahies of ^i . 

For a specific case, the values of vi and h will have been determined 
from the relations for the first stage. Let X2 be the travel at the end of the 
second stage. Then the value of /"(t'l , z) for x = X2 is given by (23), and 
22 , the corresponding value of z, can be read from the curve of Fig. 
for the value of vi applying. Then the time for the second stage, t^ , is 
given by 22(/r: + ta:) . For 2 = Z2, (20) gives !'2 , the value of w at the end of 
the second stage. This in turn gives the corresponding value of the flux 
^31^1 , and, from (21), the pull F at the end of this stage. The velocity at 
the end of this stage is given by (22), with/iCwi , 2) read from the curves 
of Fig. 6 for 2 = 22 . 

THIRD STAGE 

The total operate time is h + fi -\- h , where ^3 is the time for the third 
stage. In the relay ca.se, this is the stagger time between the operation 
of the first and last contacts. A fii'st, and frequently adequate approxi- 
mation to Is , is given by assuming the velocity in the third stage to be 
constant at the value attained at the end of the second stage. For a 
more exact determination, particularly in determining the final velocity 
for complete operation, it may be assumed that the flux in the third 
stage is constant at the value v^ifpi attained at the end of the second stage. 
The mechanical output 7 + T, in the third stage is then given by 
equation (11), with?' taken as ('2 , and .I'l and a:2 replaced, respectively, by 
X2 and X'i , the latter denoting the travel at the end of the third stage 
(zero for complete operation). Knowing the spring load 1^ between .T2 
and .Ta , h and the final velocity can be computed from the increase in 
kinetic energy T on the assumption of uniform acceleration. 
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Fig. 6 — Relations for initial motion in operation. 



FACTORS AFFECTING OPERATE PERFORMANCE 

The dominant term in the operate time is the time h for the second 
stage, which also determines the velocity and pull leyel in the third 
stage, and hence the final velocity. As U = zi{ic + i^, equation (23) 
gives the following expression for it : 
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3 ^ 2)nAiiSii{xi - .Ts) tc + is zl ^24) 

PR tc /sCwi , Zi) ' 

As can be seen in Fig. 6, j^ivi , z) is, to a rough first approximation, 
proportional to z^. To the extent that this approximation holds, ti is 
proportional to the cube root of tuAiiRiixi — X2)/{I~R), as for the last 
term in (15) and (16), except that (.ri - Xif is replaced by Ai,(Ri{xi — X2). 
From (4A) and (18), Ai6ii is given by: 

(0 — l)a;o 

Here 3^0 = A<Ro , and (2 is a minimum for that pole face area A for which 
Ai(Ri is a minimum. This condition, as determined by equating the 
derivative of the preceding expression with respect to Xq to zero, is 

given by: 

Xi = C lXq , 

or Xi/A =\/Cl-(Ro, where A is the optimum pole face area for fast 
operation. This is a smaller pole face area than that for maximum sensi- 
^tivity, for which Xi/A should equal 3to , as shown in the companion 
article^ cited above. If the pole face area satisfies the preceding condi- 
tion, the expression for Ai reduces to: 



and the term AiCRi in (24) reduces to Xi{xi — X2) multiplied by the factor 
{■\/C~L + 1)/(V'C^ " !)■ This leakage factor, which increases with the 
ratio of leakage to useful flux, is then the only term in (24) which varies 
with the design of the electromagnet. 

If fi{vi , z) were strictly proportional to 2^ equation (24) would be 
nearly independent of the coil constant N'^/R, to which tc is proportional. 
As fiivi , z) departs from this cube law relation, and also varies with vi, 
h is not independent of the coil constant. The optimum coil constant is 
that which minimizes ^1 + (3. Tn any specific case, a succession of values 
may be assumed for Gc , and h and I2 evaluated by the relations given 
above. The resulting relation between /i + (2 and Gc is similar in character 
to that described above in connection \\-ith the two stage approximation. 

6 Design for Fast Operation 

The approximations di-scussod above pro\dde a means for estimating 
the operate time attainable in specific cases, and for selecting certain 
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design variables to obtain maximum speed. In particular, it has been 
shown that for a given relay, with a specified load and power input, a 
winding can be selected to minimize the operate time. Aside from this, 
the relations show that for a given relay, the operate time depends only 
upon the load and travel and the power input, varying inversely as the 
latter at low levels of power input, and inversely as its cube root at 
high levels. 

Further conclusions can be drawm from these relations with reference 
to development studies of relays and other electromagnets. The spring 
load and travel may be considered as fixed requirements, so far as the 
magnet design is concerned. The available powder may be fixed by circuit 
considerations, or it may be related to the design by a requirement 
that the winding dissipate this power in the holding intei-val, a condition 
that imposes a minimum size on the coil and the magnet stnicture. 
Subject to this and some other limitations, there is a design choice of 
the dimensions and configuration which determine the magnetic circuit 
constants, the mass of the armature, and the eddy current conductance. 

The preceding discussion has sho\\Ti that, with an optimum choice of 
pole face area, the magnetic characteristics affect the time only with 
respect to the leakage factor, the ratio of leakage to useful flux. This 
factor may be reduced by using a "tight" magnetic circuit, but if this 
is done the factor tends to vary directly as the length of the magnetic 
path and inversely as the separation of the core and return members in 
relation to their cross sections. The leakage may be minimized by using 
a square outline for the magnetic path. The optimum speed magnet 
then has a specific configuration in which all dimensions are fixed in re- 
lation to the cross section of the magnetic members. This dimension then 
determines the mass of the armature. 

For this optimum configuration, the power, the spring load, the mass, 
and the travel determine a level of pull which minimizes the operate 
time. This pull requires a certain armature cross section if saturation is 
to be avoided, as the effective pole face area is fixed in relation to the 
cross section. If the resulting armature mass is minor compared with the 
mass of the load, the operate time attainable varies with the cube root 
of the applied powder, as in the cases discussed above. 

With increased power, however, the optimum pull and armature cross 
section increase, and must eventually reach a level where the armature 
mass becomes the dominant portion of the mass term. As this condition 
is approached, any increase in power is offset by an increase in armature 
mass, such that a lower limit is imposed on the operate time proportional 
to the travel, corresponding to an upper limit to the average armature 
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velocity attainable ^\^ith a neutral electromagnet. This upper limit is of 
the order of 100 cm/sec. Thus, for example, an operate time less than 
0.25 miUisec cannot be attained with an armature travel of 10 mil-in, 
no matter how large the steady state power applied. 

7 Release Waiting Time 

-operation-,-Telease4s-made up of an initial stage of flux change 
with the armature at rest, followed by a stage of armature motion, and 
the total time is the sum of the waiting time and the motion time. In 
release, the waiting time is usually larger than the motion time. 

There are three distinct circuit conditions under which release occurs. 
These are: 

Normal, or unprotected release, in which the coil circuit is open, and 
the only magnetomotive force maintaining the field is that of the eddy 
currents. 

Protected release, in which the coil circuit is closed through a protec- 
tive shunt, usually comprising a condenser and a resistance in series. 
The magnetomotive force comprises that of the eddy currents and that 
of the coil circuit transient. 

Slow release, in which a sleeve or short circuited winding is used to 
maintain the field and delay release. The magnetomotive force is pre- 
dominantly that of the sleeve or ^vinding current. 

The slow release case is the only one of the three for which the flux 
decay relation is accurately represented by equation (2). In protected 
release, the coil current transient is controlled by the condenser, as dis- 
cussed below. In normal release, the variation in the field linked by 
different eddy current paths results in the changes in the field pattern 
discussed in Section 1. As noted there, equation (2) applies to this 
case only as a very crude first approximation. 

In all three cases, the relation between <p and ff applying is that for 
decreasing magnetization, as illustrated by the curve for a: = in Fig. 7. 
Unlike the linear relation applying in operate, corresponding to the 
constant reluctance given by (4), the decreasing magnetization curve 
has a hyperbohc character, and is asymptotic to the saturation flux tp" . 
Residual magnetism results in a residual flux v'o , the intercept of the 
decreasing magnetization curve on the ip axis. An analytical treatment 
of the decreasing magnetization curve is given in a companion article," 
w^here it is shown to provide a satisfactory basis for predicting the release 
time of slow release relays, the third case above. The other two eases, to 
which the foUomng discussion is confined, involve both non-linear mag- 
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MAGNETOMOTIVE FORCE, 3= — »- 

Fig. 7 — Field energy relsitions in release. 

netizatioii and a variable magnetic field pattern. In an approximation 
neglecting the latter effect, no advantage is derived from an accurate 
formulation of the former, and a linear approximation to the demagneti- 
zation curve may be employed. 



NORMAL UNPROTECTED RELEASE 

As indicated above, a rough approximation to the decay of the gap 
flux, and hence of the pull, may be obtained by assuming that this flux 
decays in conformity -\\-ith equation (2), with G = Gb , where the eddy 
current conductance Ge has the same value as in operate. As a Ihiear 
approximation to the demagnetization curve, the reluctance (R can be 
taken as the closed gap reluctance 51(0). Allowance may be made for 
residual magnetism by postulating a steady state magnetomotive force 
Jo = Si(0)<po ■ Then integration of (2) gives the foUomng expression for 
the time t at which the flux has decayed to the value <p from its initial 
value tpi : 



i = Ie In 



ipi — <P( i 
ip — tpo 



(25) 
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where Ie = 47rG£/(fl(0) . Unless if_^<po is near unity, ipo may be omitted from 
this expression. If this is done, the expression may be re-written in terms 
of pull F in the form: 






(26) 



where /*' is the pull at time t, and Fi the initial pull. This expression 
follows from (25), for v-u = 0, because of the proportionahty between F 
and ip. If /*' is the spring load or tension for the operated position, and 
Fi the steady state pull there, (26) is an approximate expression for the 
release time. 



PROTECTED RELEASE 

The (commonly employed method of contact protection is to use a 
condenser-resistance shunt connected either across the coil, as in Fig. 8, 
or across the contact. Except for the steady state voltage of the con- 
denser, the circuit relations are the same for the two cases. As a first 
approximation, eddy currents may be represented by the current in a 
short circuited secondary, as in Fig. 8. The effect of such a secondary is 
determined by the value of Ie representing the ratio of its inductance 
to its resistance, aa this determines the ratio of its contribution to the 
flux to the time rate of change of the total field. 

If perfect coupling and linear magnetization are assumed to apply, 
the circuit equations for Fig. 8 can be solved and expressions obtained 
for the flux, current and condenser voltage as functions of time after 
the opening of the contact. These expressions are similar to those for 
the simple C-R-L circuit without a secondary, except that they involve 
the time constant tg as well as CR and L/R, The discharge may be either 
exponential or a damped oscillation, but while the discharge of the 
simple C-H'L circuit is always oscillatory for small values of C, the 
discharge in the cinniit of Fig. 8 is only oscillatory for an intermediate 
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Fij;. S - Coil circuit with capacitative contact protection. 
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range of capacity values. For large values of C, the discharge is damped 
by the coil resistance, and for small values of C it is damped by the cur- 
rents in the secondary. 

The latter case is essentially the condition applying for the small 
capacity values used for contact protection. Initially, the flux decay is 
opposed by the magnetomotive force resulting from both the eddy cur- 
rents and the coil current. The latter dis<;harges the condenser and 
develops a charge of opposite polarity. The initial rate of flux decay is 
therefore similar to that for exponential decay with a conductance equal 
to Gb -\- jV^/R, where R includes the external resistor as well as the coil 
resistance. The effect of the condenser charge is equivalent to a con- 
tinuously increasing coil resistance, so that the decay approaches the 
rate that would apply for exponential decay with G = Ge , and the decay 
rate is then further increased by the reverse magnetomotive force result- 
ing from the coil current reversal in the subsequent discharge of the 
condenser. 

If CR is less than ts , the effect of the coil circuit is to change the charac- 
ter of the flux versus time relation, without materially changing the time 
scale of the discharge. The initial decay is delayed, while the later 
decay is accelerated. Hence the release time is increased for a heavy load, 
corresponding to a high value of <p for release, while it may be slightly 
decreased for a light load, corresponding to a small value of (p. The utility 
of the coil-condenser circuit for contact protection results from the low 
initial rate of flux decay, which holds the induced voltage across the 
contact to a relatively low value in the initial stage of contact opening, 
when the contact separation is small. 

Except initially, the predominant magnetomotive force is that of 
the eddy currents, and this results in a changing distribution of field 
intensity, as in the case of simple release. Hence the analysis of Fig. 8, 
as described above, is not applicable quantitatively, and is therefore 
not given here in detail. Qualitatively, the relations are similar, the pro- 
tected release has a flux time relation similar in time scale to that for 
simple release, with a lower initial rate of decay, and a higher rate for 
the later stage. An illustration of this effect is included in the article by 
M. A. Logan'* cited above. No analytical treatment is available for de- 
termining the differences in release time between simple and protected 
release, at least for values of CR of the same order at ts , as in the pro- 
tection networks commonly employed. 

8 Release Motion 

It was stated in Section 7 that the waiting time in release is usually 
larger than the motion time. This, however, is not necessarily or in- 
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variably the case, and interest attaches to the conditions under which the 
motion time may become relatively large. Another important aspect of 
the release motion is the final ^'elocity attained when the armature 
strikes the backstop. This impact results hi a rebound which may, in 
the relay case, result in re-actuation of the contacts: rebound chatter. 
Rebound may be reduced by appropriate mechanical design, as dis- 
cussed in an article by E. E. Sumner,^ but its amplitude is always, other 
things being equal, proportional to the kinetic energy of the armature in 
impact on the backstop. 

Fig. 9 shows the force travel relations in release, corresponding to 
the operate relations of Fig. 4. The solid line, as in Fig. 4, represents 
the spring load, which has an operated value Fc at the point marked 2. 
This corresponds to the similarly marked pomt in Fig. 7, which shows 
the corresponding tp versus 3^ relations. The field decays along the de- 
magnetization curve for a; ^ to the point 2, where load and pull are 
in equilibrium. Further decay results in a pull less than the spring load, 
and hence in a net accelerating force producing armature motion. The 
flux continues to decrease as indicated in Fig. 7, following the path 
2-3-4. The pull follows the similarly marked curve in Fig. 9, related to 
the flux-travel path by equation (5). 

The net accelerating energy is therefore represented by the area 
marked T, lying between the pull curve and the spring load. Thus only 
this portion of the energy V stored in the spring load appears as kinetic 
energy of the armature. The remainder of V is represented by the area 




TRAVELj X — *' ^^ 

Fig. 9 — Load and pull relations in release. 
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marked W, lying below the release pull curve. This energy may be 
termed the magnetic drag: it is restored to the magnetic field and dis- 
sipated in the eddy current paths or in other circuits linking the field. 
Experimental studies have shown that the kinetic energy T of the arma- 
ture in backstop impact may vary from 20 per cent to 90 per cent of V. 
Such a wide variation in armature energy has a proportional effect on 
rebound amplitude. Within certain limitations, the following analysis 
establishes the relations controlling the ratio W/V, and also serves for 
the evaluation of the release motion time. 

EQUATIONS APPLYING 

As the motion starts after an initial period of field decay, the relations 
applying are given approximately by (2) and (3), with (R and F given by 
(4A) and (oA), and(? = Ge . The use of (4A) and (SA) is justified by the 
fact that the lower portion of the demagnetization curve, which applies 
to the later stages of flux decay, is approximately linear. The assumption 
that (2) applies, withG = Ge , is justified by the fact, discussed in Section 
1, that the initial rapid decay of the field in the outer layers of the core 
is followed by an exponential decay of the greater part (about 70 per 
cent) of the initial field. Thus if ts is written for 47rGVffl-(0) in (2), this 
expression applies in the latter stages of flux decay when the armature 
is at rest at x = 0. In this case, the value of t^ is a constant, even if the 
values of Ge and (R(0) applying differ somewhat from those applying 
in the operate case. 

Taking these equations to apply, it is desired to determine the solution 
for the initial conditions applying to the release case. This cannot be 
obtained by the approximations used for the operate case, as the sim- 
plifying condition of a nearly constant reluctance does not apply: the 
change of reluctance with x is a maximum for x ^ 0. The procedure 
that has been employed is the use of the analogue computer described 
in articles by A. A. Currie^ and E. Lakatos^ to obtain solutions for a 
restricted category of cases. The equations may be reduced to a form 
adapted to such solution as follows: 

An expression for (R/(R(0) obtained from (4A) is substituted in (2). 
Writing is for 4irGE/<^{0), there is obtained: 

Cl{1 + u) 2 , tE d/ . 
C,u ^ + 2- ^ = ^' 

where u is written for x/xo , or x/{A(Ro)- By substituting in (3) the ex- 
pression for F given by (5A) there is obtained an expi'ession for ^^ 
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On substituting this in the preceding equation, there is obtained: 



"•" 2C7. dt 



0--»)'e-™l^))-- 



To reduce this expression to dimensionless form, the following sub- 
stitutions are made; 

/ is written for -r— / Fo , where Fa is the operated load, so that / 
_ ax / 

= 1 at the start of the release motion. 
T iswrittenfor 2Cl f/ts , expressing the timei as a multiple of /e/(2Cl). 
til is written for 2mA6io/Fa. Thus /.,, is the time for travel of the mass 

through the distance ASiu for a constant accelerating force Fq . 
K is written for 2Cl tl,/il. 
The preceding equation then becomes: 



(H-«)(C. + «)(/-7v^! 



d [ ,„ . .2 1 r ^^ du 



?2. 



^^Uc.^uyu-K^-^^) = o. 



(27) 



This is the form of the release motion equation to which analogue 
computer solutions were obtained. It is a third order differential equa- 
tion giving the travel, expressed as a multiple u of ^IcRo , as a function 
of the time expressed as a multiple t of lE/i^Ci)- The boundary condi- 
tions correspond to zero initial values of travel, velocity, and accelera- 
tion, so that for r = 0, (/ = du Uh = d'u/dr = 0. 

ANALOGUE COMPUTER SOLUTION 

The analogue computer gives solutions to specific cases, and a specific 
case of (27) is defined by the \alues of Cl and /.„/'s applying, and by 
the form of the relation between /and u, defining the shape of the spring 
load. To confine the cases considered to a manageable category, they 
were limited to those for which j = 1 for all values of u: the case of a 
constant spring load. 

The constant C'l , the ratio (((iu + (JiL)/t>lo , is an inverse measure of the 
relative magnitude of the leakage field. As noted in Section 2, its value 
is usually in excess of 4. Solutions were obtained for two cases: Cl = 3 
and Cl = 5. For / = 1 , and a fixed value of Cl , the remaining param- 
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eter defining a specific case is ^.k/Ze, or the forrespouding value of K 
i'2CLtlf/(^K)- For the two selected values of Cl , solutions were obtained 
for the following values of K; 0.05, 0.1, 0.3, 1.0 and 10. 

The computer set-up for this problem is indicated in Fig. 9 of the 
article by E. Lakatos cited above.^ Solutions were obtained to (27) for 
/ = 1 and the vahies of Cl and K cited above, for the given boundary 
conditions. The solutions were obtained in the form of machine drawn 
curves gi\'ing, u \'ersus t and du/dr versus t, over the range < w < 3. 
In accordance with the character of the net accelerating force, all these 
curves were smooth and monotonic. 

For the conditions co\'ered by the solution, these results give the mo- 
tion time for a given travel directly in the form of the relation between 
T and u. It is convenient to express the values of r in terms of the cor- 
responding values of t/tAf . The results for Cl = 5 are shown in Fig. 10 
as curves of t/tM versus ii for various values of tn/tE ■ 

The other result of interest is the ratio of the magnetic drag W to 
the spring load energy V, or of (V — T)/V, where T is the kinetic 
energy at the end of a given travel: the total travel in an actual case. 
T equals m{dx/dl^'/2, or 2m{CiAaa-du/dTf/tl . For / = 1, V equals 
FoX, or Foi*-4(Ro . It follows that T/V is given by: 

T mAoto f2Ci\^ fdu 



(28) 



V 2FoU \ Ie } \di 
and therefore by : 

r _^ ^ (dju. 

V ~ 2^ V/r 



Corresponding values of u and d-ii/dr were read from the computer 
results and substituted in (28) to determine T/V and thus W/V. The 
resulting values of W/V for the case C^ = o are shown in Fig. 11 plotted 
against t^t/tR for various values of w. 

DISCUSSION OF COMPUTER RESULTS 

The significance of the results shown in Figs. 10 and 11 can be more 
readily grasped by reference to representative values of the parameters 
involved. For relay electromagnets, representative values of (51n and A 
are 0.04 cm^^ and 1 cm^, respectively, for which ASio = 0.04 cm, or 16 
mil-in. This distance is, for this case, the travel for which i( = 'i. For 
this value of Afllo , an effective mass m of 10 gm, and an operated load 
Fn of 2 X 10'' dynes (200 gm wt), t ,xf = 2 milliaec. If the load were 
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Fig. 10 — Release motion time roliitions, 

doubled and the mass halved, i^ would be halved. These vahies of tu 
arc of the same order but smaller than the values of is generally apply- 
ing, corresponding to the range in iulis covered by the figures. Similarly, 
the total travel usually lies in the range covered by the values of »., 
corresponding to travels up to 48 mils for the above value of .4(Ro . 

The curves of Fig. 11 show that the fraction ^V jV of the spring load 
energy absorbed by magneti(^ drag varies from 15 to 80 per cent over 
the range covered, and consequently that the kinetic energy at the end 
of the stroke, which determines the rebound amplitude, \'arics from 85 
per cent to 20 per cent of Y , or by a factor of 4 to 1. 

The way in which W jV varies with it and with ^w/Ze is readily luider- 
stood from physical considerations. Referring to Fig. 9, it is apparent 
that for a given time rate of field decay, the faster the rate of armature 
motion, the higher is the restraining pull and the less the kinetic energy 
T. For a short travel, or low value of », the restraining pull will be 
larger than for a high value of ii\ hence W jY decreases as xi is increased. 
As tst measures the time for a travel of .4(fto , or it ^ 1, for a force F^ , 
ImIIe is ^ measure of the motion time relative to the time of field decay. 
Hence a large value of Imi'Ie corresponds to a condition where the motion 
is slow relative to the field decay, giving little restraining pull. Thus 
IF/K decreases as /.if//s increases. 
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The c.fect of the magnetif drag on the motion time is apparent in the 
cur\-es of Fig. 10. The motion time for a given value of u increases as 
Im/Ie decreases, corresponding to increased drag. Thus a reduction in 
rebound by increased magnetic drag is, of course, accompanied by a longer 
motion time. It is of interest to note that the increasing acceleration 
results in an armature displacement which varies as the cube of the 
time, as in the relation applying in operate. 

The results shown are for the case Ct = 5. For Ci = 3, corresponding 
to higher leakage, the results are similar, but the drag ratio W/V values 
are from 5 to 10 per cent lower than those for CV = 5, and the values 
of t/tM are correspondingly smaller. Thus the value of Ct applying has 
only a secondary effect on the release motion. 

While these results are formally limited to the special case of a con- 
stant retractile force, the conclusions are of broader application. A load 
curve that decreases as the gap opens modifies the solution for the same 
operated load Fq only in reducing the net accelerating force, increasing 
the drag ratio and the motion time. The residual magnetism, neglected 
in this discussion, has a similar, and minor, effect. The results show Cl 
to have little effect on the motion, which is governed primarily by the 
values of u and iju/ij applying. 
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Fig. 11 — Energy absorbed by magnetic drag in release motion. 
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As minimizing rebound is of major importance in relay design, the 
relations controlling the kinetic energy to which it is proportional are of 
design interest. Most of the quantities appearing in the above relations, 
such as the spring load and the travel, are fixed by operating require- 
ments or other design considerations. There is, however, considerable 
freedom in the design relations fixing the effective mass of the armature, 
particularly in selecting its axis of rotation. The preceding relations 
show that the magnetic drag ratio can be increased, and rebound re- 
duced, by the decrease in Im resulthig from a decrease in effective mass. 
While such a decrease increases the ratio i/tn , it results in a net decrease 
in the value of t, the motion time. A low effective armature mass is 
therefore advantageous both for fast release and for a high drag ratio 
for the reduction of rebound. 
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